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We improve on Cutkosky's cutting rules which is used to calculate the contribution of 
the singularities of Feynman propagators to Feynman amplitude. The correctness of the 
improved cutting rules is verified by the calculations of the conventional loop momentum 
integral algorithm. We also find that only the process of the intermediate on shell particles 
that is same as real physical process can exist in quantum vacuum. 
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I. INTRODUCTION 

Cutkosky's cutting rules is a simple algorithm to calculate the contribution of the singularities 
of Feynman propagators to Feynman amplitude [1], and has been widely used in calculating the 
imaginary parts of physical amplitudes. Many relevant problems, for example the Abramovskii- 
Kancheli-Oribov cu..i„. Rule. the cutti„. in finite te.pe.a.n., field theory 1,1. the 

unitarity method in the cutting rules [J], the Feynman Tree Theorem [5!], the dispersion method 
to calculate Fe.a.aa a.pl.tudc. fl. the du^ity tclat.on to calculate Feyn.an integral. the 
cutting rules and the unitarity of noncommutative field theory from string theory [8] , and the Hopf 
algebra in the cutting rules [sl, have been widely discussed around Cutkosky's cutting rules. From 
these discussions one finds that Cutkosky's cutting rules is a very useful method both in actual 
calculations and in searching for underlying physical meanings of Feynman diagrams. 

But all of the discussions (including Cutkosky's initial discussion) haven't carefully investigated 
the algorithm's structure, so haven't noticed the immaturity of the algorithm and the relevant 
underlying physical meanings. Here we give a careful investigation about it. Consider the Feynman 
propagator, 

Df{x -y)= I e-^P<--y^ , (1) 



there are two singularities po = ±(-\/p^ + m"^ — is) in the denominator. Since the cutting rules is 
proposed for calculating the contribution of these singularities to Feynman amplitude, it should 
point out that in what cases which or both of the singularities have contribution to Feynman 
amplitude, since in some cases only one of the singularities has contribution to Feynman amplitude. 
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for example in the one loop self energy, 



(27r)^ k? — m\ + ie {k — pY — ra^ + ie 



where the outline momentum p is on mass shell pq = + m? (m is the mass of the outline 
particle), there are four singularities in the propagators, if keeping all of the contributions of the 
singularities, we have the result by Cutkosky's cutting rules: 

r A./ ^ T f d^k -i 1 

ImM{p^p) = Im / -TTT-UTo 2~r~Ti ^2 TT~ 

J (27r)* k^ — mf + le [k — pY — ra2 + ie 

(m^ + m\ + ml — 2m?m'i — 2m?m'2 — 2rn\'m'2)'^l'^ 

X [9[mi — m — 7TI2] + 0[m — mi — m2\ + 9[m2 — m — mi]) , (2) 

where 9 is the Heaviside step function, there are three Heaviside step functions (there should be four 
step functions, but one of them has been eliminated by the energy conservation law in each vertex), 
obviously such result is wrong, since according to the Breit-Wigner formula ImM{p — >■ p) = mV (F 
is the decay width of the outline particle) [10;] only the second step function has contribution to the 
le self energy. There are some discussions about the problem of the selection 
^, y, ^, I9I, but a thorough and practical prescription hasn't been present. 
So setting up an explicit and practicable algorithm for finding the singularities having contri- 
bution to Feynman amplitude becomes very important. In order to solve this problem we firstly 
discuss the origin of the three step functions of Eq.(2). It is obvious from the calculation of Eq.(2) 
that the three step functions come from the three different combinations of the singularities, as 
shown in Fig.l. 

mi 

m / 
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FIG. 1: The three different combinations of the singularities of the one loop self energy diagram. 



In Fig.l we stipulate that an arrow in a propagator denotes the propagator is on mass shell, and the 
arrow direction denotes the particle of the propagator incoming into or outgoing from the vertex of 
the Feynman diagram as convention, and the energy component of the momentum along the arrow 
direction is positive. Thus an arrow in a propagator represents a singularity of the propagator, and 
the direction of the arrow tells us which singularity it represents. Since we call a kind of combination 
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of propagators' singularities as a 'cut', the three diagrams in Fig.l represent three different cuts of 
the one loop self energy diagram. Prom the calculation of Eq.(2) we find that the first, the second 
and the third cuts of Fig.l contributes the first, the second and the third step functions of Eq.(2) 
respectively. From the previous discussions we know only the second cut of Fig.l has contribution 
to the one loop self energy amplitude, the other two cuts of Fig.l must be eliminated. Comparing 
the three cuts of Fig.l we find the difference between the second and the other cuts is whether the 
arrow directions are along a contrary circumrotating direction in the loop: in the second cut the 
arrow directions are along a contrary circumrotating direction in the loop, but in the other two 
cuts the cases are just opposite. Since in physical meanings the arrow direction stipulated here 
represents the propagating direction of the intermediate on shell particle in quantum vacuum, the 
arrow directions being contrary in a loop means the corresponding intermediate on shell particles 
propagating from a same vertex to another same vertex. If we can generalize this conclusion to 
the Feynman diagrams with arbitrary number loops and arbitrary number outlines, we will draw 
the conclusion after some simple analysis that all of the intermediate on shell particles in quantum 
vacuum propagate along a same direction: from the incoming particles to the outgoing particles of 
the Feynman diagram. 



We assume that the above-mentioned generalized conclusion is right, i.e. only the cut in which all 
of the intermediate on shell particles propagate along a same direction-from the incoming particles 
to the outgoing particles of the Feynman diagram-has contribution to the Feynman amplitude. 
This means only the cut in which the arrow directions stipulated here of the cut propagators (i.e. 
the propagating directions of the intermediate on shell particles) cannot form a closed circle in any 
of the loops has contribution to Feynman amplitude. 

In order to suitable for actual calculations we need to solve a potential problem in a kind of 
Feynman diagrams part of whose innerlines are inserted by self-energy or tadpole diagram. This 
kind of diagrams, for example see the diagram of Fig.2 



II. IMPROVE THE CUTTING RULES 




FIG. 2: A diagram one of whose innerlines is inserted by self-energy diagram. 
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In Fig. 2 if A and B propagators can be cut simultaneously, it will contribute an infinite large term 
to the Feynman amplitude when A and B particles have same masses, since there will be two 
same delta functions in the loop momentum integral. So such cut should be prohibited. In exact 
words, beyond one cut among the propagators which are connected each other by self-energy or 
tadpole diagrams is prohibited. Under this constraint the contribution of the cuts of Fig. 2 will be 
equal to f{in\)/{m\ — m^) + f{m'g)/{m\ — m\), where / denotes the amplitude of cutting A 
propagator and removing the denominator of B propagator or the contrary case, and mA-,fnB is 
the mass of A and B propagators respectively. When rriA = niB = m, the cut contribution will be 
equal to df{m?)/drn?. Through some simple deduction we can draw the conclusion that all of the 
contributions of the cut diagrams which have arbitrary number propagators which are connected 
each other by self-energy or tadpole diagrams are finite under the above-mentioned constraint. To 
specific, if the number of the propagators connected each other by self-energy or tadpole diagrams 
is n and all of the masses of the connected propagators are same, the contribution of the cut 
diagrams to the Feynman amplitude will be equal to d^h{m?) / (dm?)^ , where m is the same mass 
of the connected propagators and h is the amplitude of cutting one connected propagator and 
removing the denominators of all of the other connected propagators. 

Now we need to determine the contribution of the singularity of a propagator to Feynman 
amplitude. Base on the mathematic formula 

^ P-^iTT5{a), (3) 



a±ie a 

where P denotes the Cauchy principle value, the contribution of a determinate singularity of a 
propagator, i.e. a cut of the propagator, to Feynman amplitude should be 

3^-— ^ -i7re\po] <5(/ - m^) , (4) 

where the propagating direction of the cut propagator has been given, and the term 9[po] is just 
used to guarantee the energy component of the momentum p along the given propagating direction 
is positive, as required by the previous stipulation. Since the delta function has two singularities 



o{p -m) = + 



PO + \/p^~+m2 Pq - + 771,2 

we find that Eq.(4) is different from Cutkosky's formula: there is only one singularity in Eq.(4), 
but in Cutkosky's formula both of the singularities exist [ll]. We can simplify the expression of 
Eq.(4) as follow: 

y -iTr6{po - Jv^ + m?)/ {2 J p2 + m?) . 

pz _ _|_ V V 
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Furthermore we notice that when we integrate a loop momentum, all of the selected singularities 
in the loop will be calculated i times where i is the number of the selected singularities in the loop. 
This is because there are i times integrations for the i number delta functions, so all of the selected 
singularities will be calculated i times. Therefore if a loop has i number cut propagators, the result 
should be multiplied by an extra factor i. For a complex Feynman diagram the number of loops 
can be large than one and the loops can be entangled each other, so in order to avoid repeated 
counting we only need to count the extra factor of the loops which don't contain any smaller loop 
in it. 

In the view of mathematica the principle value and the single delta function in Eq.(3) has 
similar contribution to Feynman integrals, both contribute normal results, only two or beyond two 
propagators' delta functions in a momentum loop have abnormal contribution: the discontinuity 
of the Feynman amplitude. Thus a Feynman amplitude can be divided into three parts: the 
contribution that all of the Feynman propagators are replaced by their Cauchy principle values, the 
contribution that part of the Feynman propagators are cut (i.e. replaced by the delta functions) but 
none of the momentum loops are cut beyond once, and the contribution that there are momentum 
loops which are cut beyond once. Obviously only the last part contributes the discontinuity to 
the Feynman amplitude. Since the purpose of proposing the cutting rules is to calculate the 
discontinuity of any Feynman diagram |l|, the algorithm should reach such goal: all of the 
discontinuity part of a Feynman diagram are denoted by cut, and the remaining part aren't denoted 
by cut. According to such requirement the cutting rules should satisfy the following two conditions: 
1) the least number of the cut propagators in a cut momentum loop is two; 2) the uncut propagators 
in a cut momentum loop are replaced by their Cauchy principle values, the propagators in an 
uncut momentum loop keep unchanged, but their contribution to the discontinuity of the Feynman 
diagram needs to be removed. 

Summing up all of the above-mentioned discussions we get the following improved cutting rules: 

1. Cut through the Feynman diagram in all possible ways such that the cut propagators can 
simultaneously be put on mass shell and the cut propagators' propagating directions (it is 
required that the energy component of the momentum along the propagating direction of the 
cut propagator is positive) cannot form a closed circle in any loop; keep the number of the cut 
propagators in a same propagator chain (where the propagators in a same propagator chain 
means all of the propagators are connected each other by self-energy or tadpole diagrams) 
less than two; keep the least number of the cut propagators in any cut loop is two. 
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2. For each cut propagator with definite propagating direction, replace — m? + ie) — >■ 
—m5{pQ — + m?)/ (2-\/ p-^ + m^) where the momentum p = (pO)P) is along the given 
propagating direction of the cut propagator; for each uncut propagator in a cut loop, apply 
Cauchy principle value to it; for the propagators in an uncut loop, keep them unchanged, 
but their contribution to the discontinuity of the Feynman diagram needs to be removed; 
for each cut loop which doesn't contain any smaller loop in it, multiply by factor i where i 
is the number of the cut propagators in the loop, then perform the loop integrals. 

3. Sum the contributions of all possible cuts. 

We give an example to illustrate the improved cutting rules. Fig. 3 shows the cuts contributing 
real part to the two loop self energy according to the improved cutting rules. From Fig. 3 we can 
see several different processes happen in the quantum vacuum that the incoming particle decays 
into some intermediate on shell particles then these intermediate on shell particles decay into the 
outgoing particle. In addition we find that the singularities of Feynman propagators also contribute 
real part to Feynman amplitude. 






FIG. 3: Cuts contributing real part to the two loop self energy. 



III. CORRECTNESS OF THE IMPROVED CUTTING RULES 

In this section we will examine the correctness of the improved cutting rules. We will calculate 
the imaginary parts of some Feynman diagrams by two methods. One method is the improved cut- 
ting rules, the other method is the conventional integral algorithm, i.e. the Feynman parametriza- 



12l |. We will compare the two methods' 



tion, wick rotation and the dimensional regularization 
results to judge whether the improved cutting rules is correct. In the following calculations we 
have used the program packages FeynArts and FeynCalc [l^. The examples of calculating the 



contributions of the singularities to real parts of Feynman amplitudes can be found in Ref. 



l4. 
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Firstly we calculate the imaginary part of the two-loop two-point Feynman diagram shown in 
Fig.4, where the cut has been drawn and = m?. 

mi 




FIG. 4: Cut contributing imaginary part of a two-loop two-point diagram. 



Using the improved cutting rules we obtain 

d'^ki d^k2 1 1 



Im 



(27r)4 {2-kY kl-ml + ie kl-ml + ie {ki + k2 - v? -mj + ie 
= -^J d^hd% e[kio] Sikj - ml) e[k2o] 5{kl - ml) 0[po - ^lo - k2o] 5{{p -ki- A^a)' - ml) 



m —3m-, 



1 / ^ dx 
64m7r^ Jmi 



111 
m\\m — m 



2\2 



V2, 



(mx -\- m?)^ T, — — ' Q\m — Smil . (5) 

^ ^' m^ + ml-2mxi i w 



Using the conventional integral algorithm we also obtain 
d^ki d^k2 1 1 



Im 



{27r)^ (2it)^ kl-ml + ie kl-ml + ie {ki + k2 - p)'^ -ml + ie 

1 fi /"i 2 

= jIm dx (i2/(A-2m^y)(--27-|-lnl67r^-|-lnj/-lna;(l-a;)-lnB)lnB 

2567r* Jo Jo e 

= -—^1 dx I dy {A - 2m^y){- - 2-i + h).lQT:'^ + \nij-\nx{l-x)-2\D.\B\)e\-B] 
2567r"^ Jo Jo e 

= 5- / dx{— ^ — m?ln — 6* m - 3mi , (6) 

where e = 4 — D, 7 is the Euler constant, and 



A 2,2 "^1 

A = m +m^ , , 

x{l — X) 

B = m'^y'^ — Ay + ml , 

m — mi ^ (m^ — 2mmi — 3m? )^/^ 

^1,2 = ^7 X • (7) 

2{m — mi) 

In Eq.(6) we have used the formula Ina = In |a| — ztt 0[— a]. Both Eq.(5) and Eq.(6) are too complex 
to be expressed by analytical formula, so we only compare their numerical results as shown in Fig.5. 
From Fig.5 we find that the results obtained by the two methods coincide very well. In order to 
show to what extent the coincidence is we calculate the relative discrepancy of the two results, as 
shown in Fig.6. We note that we have used the software Mathematica 5.0 to calculate the result, 
and we have set the precision goal as 10~^°. The maximum relative discrepancy of Eq.(5) and 
Eq.(6) is 3.4 x 10^^^ which is less than the set precision goal, so we can draw the conclusion that 
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ImAI (m -> m) 

m2 




0.05 0.1 0.15 0.2 0.25 0.3 m 



FIG. 5: Numerical results of Eq.(5) (denoted by dots) and Eq.(6) (denoted by line). 
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FIG. 6: The relative discrepancy of Eq.(5) and Eq.(6). The amplitude with a subscript "cut" denotes it is 
obtained by the improved cutting rules, and the amplitude without subscript is obtained by the conventional 
integral algorithm. 

Eq.(5) is equal to Eq.(6) at the precision 10"^'^. Besides, we haven't shown the data whose absolute 
values are greater than 10~^^ in Fig. 6 in order to make the great majority of the data clearly shown 
in the graph. 

Then we calculate the imaginary part of the one-loop three-point Feynman diagram shown in 
Fig.7, where = and Pi=P2 = ml- 




FIG. 7: Cuts contributing imaginary part of a one-loop three-point diagram. 



Using the improved cutting rules we obtain 

d^k 1 



(27r)^ k"^ —m\ + ie {k — pi)^ — m\ + ie {k — pi — ^2)^ — m\ + i£ 
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P 



+ 



+ 



Stt^ [ {k—piY — m\ + ie 

e[ko]5{e - mf )g[pio - k^]5{{pi - kf - mj) 

{k — pi — P2Y — + i£ 
0[ko - pm]5{{k - pif - mf )6i[pio + P20 - ^o]'^((pi +P2- kf - mj) 



k^ 



l£ 



9[ma — 2mi 



In 



+ 



IGnmay m'^ — Amf 
6[mb — 2mi] 



2ml 



\J{ml - ^rnl){ml - 4mf) 



— 2m^ + \J (m^ — 4m^)(m^ — 4mf) 



In 



TTiamb + (m^ — 4ml) {ml ~ 4mf ) 



Svrma Y — 4m^ 
Using the conventional integral algorithm we also obtain 

r ^ ^ . d^k 1 1 

/m(— i) 



(27r)'^ fc2 — mf + ie {k — pi)"^ — m\ + is {k — pi — p2)^ — mj + ie 



1 fl—x 

^ Im / d x / 

167r"' Jo "'0 



dy 



+ ml + x y — x — y 



1 /"^ nl—x 

/ dx dy 6{mlx'^ + mly'^ + mlxy — mix — ml y + mj) 

Idtt Jo Jo 



1 



:ln 



167rmawm^ — Ami 



m" 



mf: 



2mi 



^ {ml - 4ml){ml - 4mj) 
2ml + V ("^a ~ 4m^)(m^ — 4mf ) 



mamb + (m^ — 4ml){ml — 4mj) 



.rriamb 



ml — 4ml){ml — 4mj) 



0[ 



mb 



mi] + 



167rmaJ ml — 4ml 



X In 



— 2mf^ 



m^ 



[ml-4ml){ml-4ml) jria , 

d[ mi\ U[mi J . 



2 -2ml + J {ml - 4ml){ml - 4mf) 



(8) 



(9) 



In Eq.(9) we have used the formula Eq.(3). Obviously Eq.(8) is equal to Eq.(9). Eq.(8) also 
coincides with Eq.(3.24) of Ref.Q. 

Thirdly we calculate the imaginary part of the one-loop four-point Feynman diagram shown in 
Fig.8, where r = {2p,pez) and q = {2p, —pe^). 



m 



m 



+ 



+ 



m 



+ 



+ 



FIG. 8: Cuts contributing imaginary part of a one-loop four-point diagram. 
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Using the improved cutting rules we obtain 



[ d'^k 1 



(27r)^ k"^ + ie {k — q)'^ — m? + ie k'^ — m? + ie (k + r)^ — m? + ie 



87r2- 



P / (fk 



r 6'[fco](5(fc2 - m2)g[go + rp - kQ]5{{q + r - kf - m^) 
(/c — r)'^{{k — r)2 — m^) 



e[ko]5{k^ - m'^)0[ro - ko]6{{r - kf - m^) 9[ko]6{k '^ - m'^)e[qo - ko]5{{q - kf - m^) 

~t~ I O / / 7 '. To ~t~ 



k^({k + — 



g[fco]^(fc2)g[ro - koMr - kf - rn") e[kQ]6{k^)e[qo - koWq - kf - m^) 

"T /lO n\//T \n n\ "T 



{kP' — ■w?){{k + rf — m?) 



e[p- /»/v^] 0[p - m/2] V - m'^p + 2;-;?%^ Ip- - irfi 

0[p - 2m/^l , 15p2 - + 4p^/V - I2m2 

In ^^^^^^^^^= 

1287rm2p2 ]^5p2 _ _ 4p^Qpi _ I2m2 

Using the conventional integral algorithm and Eq.(3) we also obtain 

ct^k I I 1 1 



+ 



(10) 



1 



167r2 



(27 

1 



+ ie {k — qf — rm? + ie k'^ — m? +ie {k + rf — m? + ie 







l-x 



dy 



l~x-y 



dz 



J rl pl—x 

= ifi 2 2 -^"^ dx dy 
ibiT^m^ Jo Jo 



(3p2 x2 + 3p2 y2 _ ]^gp2 xy — 3p2 x — 3p2 y — m? z + — ie) 



3p2 a;2 + 3p2 y2 _ jQp2 xy -\- (jn^ — 3p2)a; + [vrP- — 3p2)y — ■je 



3p2 2;2 _|_ 3p2 y2 _ ]^gp2 x y — 3p2 x — 3p2 y + — ie _ 

= Jo ^^Jo ^'^ l*^^^^^ ^ ~ "^^^^ X y + (m^ - 3p^)x + (m^ - 3p2)y) 

- 5(3p^ + 3p2 - 10p2 a; y - 3p2 a; - 3p^ y + m^)] 

— In 



1287rm2p2 g^s _ ^2^ _ 2'm?\/Ap'^ — 

rn ]Q\2ni .,, „ , o , o — y 



H T^ri: ;t-t; in 



+ 



1287rm2p2 "' 9(9^8 _ ^2^ _ 2m'^^/4p^ - vn?) 
- ^] {^p^ - m^p + 2m2V4p2-m2)(15p2 + 4pV V - ISm^ - 4m2) 



(11) 



1287rm2p2 g^g^s _ _ 2w?^/Ap'^ - m?){lhp'^ - Ap^/Qp^ - I2m? - Am?) ' 
Obviously Eq.(lO) is equal to Eq.(ll). 

Lastly we calculate the imaginary part of the two-loop two-point Feynman diagram shown in 
Fig.9, where p^ = M^. 




FIG. 9: Cuts contributing imaginary part of a two-loop two-point diagram. 
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Using the improved cutting rules we obtain 

d^ki d^k2 1 1 



J (27r)^ {2tt)D kf + iekl-m? + ie 
1 1 1 



(^2 — kiY — -\- ie kf — w? + ie {k\ — p)^ — m? + ie 
= T 2epfd^d^r 1 

""^ J (27r)^(27r)«Kp-A:i)2((p-A:i)2-m2) 
X 4iTT^ e[kio]5{kl - m^)9[k2o]5{k'l - m^)e\po - kw - k2Q]6{{p - ki - fcs)^ - m^) 
_ 27r^e[kw]5{kl)e[po - k^oMp - fci)2 - m^) 
{k2 — m?){{k2 — kiY — w?){k'l — w?) 

2TTH[kw]5{kl - m2)%o - kio]5{{p - k^f - w?) 



kl{kl-w?){{k2-kiY -w?) 

roo poo rl 

—r^P / kldki / kldk2 / dx 
647r-^ Jo JO J-i 

6{M - ^Jm? + kl - ^Jrln? + kl - ^Jw? + kl + kl + 2kik2 x) 

w? + k\ + k2 + 2kik2 X 

1 



(M2 + m2 - 2M^m2 + A;f)(M2 - 2M^m^ + kj) 

2^2^2ep(g/2) f d^ki i.,,2x/ir ,1../ , ^2 2 

+ (J)D/2^2 J J^G[kio]S{kl)e\po - k^oMp - k,f - 
2^2^2ep(g/2) /•! dx f d^ki 



(47r)^/2m2 Jo (x2 - a; + 1)^/2 7 (27r) 

0[M-3m] 2 f^^^^ v/(x2 - 1)(M2 - 2mMx - 3m2) , 

— Tft I dx 

647r3M J I (M2 + m2 - 2mMx)3/2(M - 2mx) 

e\M-m] ,,,2 2n/a M2-m^ 

— , „ ^ VM2 -4m2(A-ln ^ h4 ^ , (12) 

2567r3m2M ^ m2 ' ^ ' 

where 

2 m2 
A = - + 21n47r-27-21n^. (13) 
e 

Using the conventional integral algorithm and Eq.(3) and the formula Ina = ln|a| — md\—a\ we 
also obtain 



Im ij?^ 



d^ki d^k2 1 1 

(27r)^ (27r)^ fc2 + _ ^2 + 
1 1 1 



(^2 — fei)^ — w-^ + ie kf — m? + ie (ki — p)^ — m? + ie 
= 2^^'"^ L '''' L L Wy 



dz r . , 1 — v — 2; 

A + In ■ 



/2 + {w? — M^)y + m22; — 



a; — a;2 



12 



- H—r2 — + y + ^ + -r— if 



5 ze) - ln( ^ ^ + y + 2: - 



^/ / dy dzS(M'^y'^ + (rn^ - M'^)y + rn'^z) A + ln 

256n-^ Jo Jo Jo L 

, , M2y2 M^y 1-y- z.i 
In ;5 TT + y + z-\ ^ 



1 - y - ^ 

X — x"^ 



+ 
+ 



'Gtt^ Jo ^ Jo ^ Jo ^ M'^y' 



l-y mM^_MV 1=^] 

J L 771/ rm. ^ x—x ^ 



2567r3 

1 



-Im I dx dy 



Jo 



Jo 



P- + (jv? — M'^)y + m?z 

M^y^ M^y ^M'^y^ M^y . ■ 
K TT + y - le) - In { 7^ ^ + l — ie) 



= 9A^9 ^ M^-m^) A- 21n ti— + ^) 



zooTT'^in'-M^ mm 
+ ^[^-fft^r;-^] Um^ - - Mv/M^^)A + 2(M^ - m^)(l - In 

, ,0 ,0 



+ 



g[M - 3m] 
^ 2567r3m2M2 



M 



/Lf2 47772 „ 

• - MV'M2 - 4to2)A + MV'M2 -4m2(ln + - 4) 



(14) 



where 



F = 



M^y^ - M^y - ^mPy + 4m2 \ V2 



(15) 



, y(M2y-M2 + m2; 
M2 + ^ - 10m2M2 + Qm^ 
^^'^ ~ 2M2 

Both Eq.(12) and Eq.(14) are too complex to be expressed by analytical formula, so we only 
compare their numerical results, as shown in Fig. 10. From Fig. 10 we find that the results obtained 
by the two methods coincide very well. In order to show to what extent the coincidence is we 

ImM (M^M) 



0.00015 
0.0001 

. 00005 



-o.oooosf 



FIG. 10: Numerical results of Eq.(12) (denoted by dot) and Eq.(14) (denoted by line) with A = 0. 



calculate the relative discrepancy of the two results, as shown in Fig.ll. Similarly as Fig.6 we have 
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ImMcut (M ^ M) - ImM (M ^ M) 
ImM (M ^ M) 

4 • 10"^^ 




FIG. 11: The relative discrepancy of Eq.(12) and Eq.(14). The amplitude with a subscript "cut" denotes it is 
obtained by the improved cutting rules, and the amplitude without subscript is obtained by the conventional 
integral algorithm. 



used the software Mathematica 5.0 to calculate the result, but we have set the precision goal as 
10"i2_ As we see from Fig. 11, the maximum relative discrepancy of Eq.(12) and Eq.(14) is less 
than 10^^^ which is less than the set precision goal, so we can draw the conclusion that Eq.(12) is 
equal to Eq.(14) at the precision 10^^^. Besides, we note that a similar discussion has been present 
in Ref.Q. 

There are a lot of similar discussions which are less typical, so we don't list them here. All 
of the discussions show that the results of the improved cutting rules coincide with those of the 
conventional integral algorithm very well. So we think that the improved cutting rules is correct. 



IV. DYNAMIC LAW OF THE INTERMEDIATE ON SHELL PARTICLES IN 

QUANTUM VACUUM 

We have mentioned in section two that the improved cutting rules is based on the assumption 
that only the cut in which the intermediate on shell particles propagate along a same direction-from 
the incoming particles to the outgoing particles-has contribution to Feynman amplitude. What 
can we get from the assumption? It is obvious that this assumption guarantees two things: 1) 
the intermediate on shell particles cannot create themselves circularly in quantum vacuum, for 
example the cases of the first and the third diagrams of Fig.l are forbidden, thus the energy of any 
intermediate on shell particles is less than the total energy of the incoming particles; 2) the time 
order of all the vertices connected by the intermediate on shell particles is definite, since along the 
propagating direction of the intermediate on shell particles the time increases and the intermediate 
on shell particles propagate along a same direction-from the incoming particles to the outgoing 
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particles (for example the time order of the vertices in Fig. 3 is definite, and the case of indefinite 
time order of the vertices connected by intermediate on shell particles, like the cases of the first and 
the third diagrams of Fig.l, is forbidden), thus the time order of the processes of the intermediate 
on shell particles in quantum vacuum is definite, i.e. beginning from the incoming particles lastly 
ending at the outgoing particles. 

On the other hand, from all of the cutting diagrams of this paper we find that a correct cutting 
must cut a Feynman diagram into several completely separated parts (this can be guaranteed 
by the constraint in the improved cutting rules that the least number of the cut propagators 
in any cut loop is two). Combining all of the above discussions we can get such a picture in 
mind about processes of the intermediate on shell particles in quantum vacuum: the incoming 
particles decay into intermediate particles, then the intermediate particles decay into the outgoing 
particles; or the incoming particles decay into intermediate particles and part of the outgoing 
particles, then the intermediate particles decay into the remaining outgoing particles; or part of 
the incoming particles decay into intermediate particles, then the intermediate particles and the 
remaining incoming particles decay into the outgoing particles; or the other intermediate processes 
at two level and beyond; all of the creation, propagation and annihilation of the intermediate 
particles are along a same direction: from the incoming particles to the outgoing particles, while 
the circular creation of the intermediate particles is forbidden. Thus we find that the processes of 
the intermediate on shell particles in quantum vacuum are same as real physical processes, the only 
difference between them is that the final states detected by experiments are the outgoing particles 
of the real physical processes, not the intermediate on shell particles, so we haven't noticed the 
existence of the processes of the intermediate on shell particles in quantum vacuum. 

In fact in some special physical processes the process of the intermediate on shell particles can 
be observed. For example, when a physical colliding process happens through the way that the 
incoming particles fuse into an intermediate particle, then the intermediate particle decays into the 
outgoing particles, the curve of the cross section will present a peak at the position where the total 
energy of the incoming particles is equal to the mass of the intermediate particle [10] . Obviously 
this peak in the curve of the cross section indicates the presence of the process of the intermediate 
on shell particle in quantum vacuum. 
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V. CONCLUSION 

In this manuscript we discuss how to improve Cutkosky's cutting rules to make it more exphcit 
and more practicable. We point out that a propagator has two singularities, so the cutting rules 
must tell us which singularity is used in the calculations. In order to solve this problem we introduce 
a concept: the propagating direction of the on shell propagator, and we stipulate that the energy 
component of the momentum along the propagating direction of the on shell propagator must be 
positive. Thus the propagating direction of the on shell propagator represents the real propagating 
direction of the intermediate on shell particle in quantum vacuum. Under such stipulation which 
singularity being used will be given beforehand in a cut propagator with definite propagating 
direction. Through some careful analysis we acquire the improved cutting rules. In order to verify 
the correctness of the improved cutting rules we calculate the imaginary parts of several Feynman 
diagrams by two methods: the improved cutting rules and the conventional integral algorithm. 
We find that the results obtained by the two methods coincide very well. So the improved cutting 
rules should be correct. 

Through the foundation process of the improved cutting rules we find that the processes of the 
intermediate on shell particles in quantum vacuum must obey the following rules: 1) they separate 
a Feynman diagram into several completely separated parts; 2) all of the creation, propagation and 
annihilation of the intermediate on shell particles are along a same direction-from the incoming 
particles to the outgoing particles; 3) the circular creation of the intermediate on shell particles 
(for example see the cases of the first and the third diagrams of Fig.l) is forbidden, and the energy 
of any intermediate on shell particles is less than the total energy of the incoming particles. Thus 
the processes of the intermediate on shell particles in quantum vacuum are same as real physical 
processes, the only difference between them is that the final states detected by experiments are the 
outgoing particles of the real physical processes, not the intermediate on shell particles. In other 
words, only the process of the intermediate on shell particles that is same as real physical process 
can exist in quantum vacuum. 

The present conclusion is only a pilot study. The problem of the intermediate on shell particles 
in quantum vacuum needs more investigations. We can expect that through more investigations 
people can reveal more and more secrets of the quantum field vacuum. 
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